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Cs) , Abstract 

A systematic study of the strongly interacting matter under extreme conditions via 
the form of the thermal ratio of the disorder deviation is presented. The evolution 
of fermi- and bose- particles (quarks and gluons) is studied in the framework 
of multi-particle correlation and distribution functions to predict the size of the 
finite-temperature deconfmed phase. 
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O 1 Introduction 

^. 

Q-i! The existence of a deconfmed phase of gluons and quarks has been predicted by quan- 

Q_i! turn chromodynamics (QCD) [1]. Aspects of QCD at finite density are very important as 

well as instructive in the sense of the observation of quark-gluon plasma (QGP), especially 
from the heavy vector meson (HVM) suppression, strangeness enhancement, hadron yield 
distributions versus temperature and dilepton excess. There is large amount of consis- 
tency among the different signatures of deconfmed phase transition to a new excited state 
of matter. In the case of HVM suppression (e.g., cc bound state suppression), the CERN- 
SPS experiments clearly show that it's visible in Pb-Pb collisions starting from a definite 
value of the energy density in the right range. Enhanced strange particle production in 
Pb-Pb collisions offers another indication for the phase transition: in the deconfmed phase 
strange quarks are more easily produced. The threshold for the producing strange quarks 
is much lower than that for strange hadrons. In addition, the mass of the strange quark 
goes down in the case of the restoration of the chiral symmetry. Hence, the observation 
of large ratios for strange particles is considered as a remnant of the unconfined phase. 
The dilepton abundance, while well explained by ordinary soucers in all other cases, for 
Pb-Pb collisions a neat excess is instead observed. 

There is a very popular point of view in literature that a deconfinement phase tran- 
sition is predicted to occur at the typical energy scale involved, the temperature T c ~ 
QCD scale Aqcd ~ the pion mass ~ the mass of a strange quark m s ~ 140-200 
MeV. Putting all indications above mentioned together, the case of a high density state 
of matter onsetting at a critical temperature around 170 MeV is reasonably suitable. 

Among the issues related to QGP, we attract attention to the problem of deconfmed 
phase through the calculation of correlation and distribution functions [2] in the thermal 
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theory of quantized fields. We consider the semiphenomenological model for the decon- 
finement existence within the framework of the Langevine-type equation. To do this, we 
follow the standard theory of quantized fields replacing: 

1. the asymptotic field operators and 

2. the vacuum expectation values 
by 

1. the thermal field operators and 

2. the thermal statistical averages, 

respectively, in order to formulate correlation and distribution functions of produced par- 
ticles. 

The method of Langevin equation and its extensions to the quantal case have been 
suggested and considered in papers [3-5] and [6-8], respectively. We propose that real 
physical processes to happen in the finite-temperature QCD should be replaced by a one- 
constituent (e.g., gluon or quark) propagation provided by a special kernel operator (in 
the master evolution equation) to be considered as an input of the model and disturbed 
by the random force F. We assume F to be the external source proposed as both a 
c-number function and an operator. The master equation is an operator one, so that 
there appear new additional issues about the commutation relations and the ordering of 
operators, which do not exist in the classical case [8]. 

Based on the thermal operator-field technique, we introduce a thermal ratio of the 
disorder deviation (TRDD), reflecting the degree of deviation, from unity, of the ratio 
of the two-particle thermal momentum-dependent distribution to two one-particle ther- 
mal distribution functions of produced particles, gluons and/or quarks (g/q) in a partly 
deconfined phase state. We study the four-momentum correlations of identical particles 
which can be both useful and instructive to infer the shape of the particle emitter-source. 
We estimate the sensitivity of the TRDD-functions to the size of the emitter. Within 
these features, the canonical formalism in a stationary state in the thermal equilibrium 
(SSTE) is formulated, and a closed structural resemblance between the SSTE and stan- 
dard quantum field theory is revealed. 

To clarify the internal structure of the disordering of particles, we use the consistent ap- 
proach based on the evolution of dynamical variables as well as the extension to different 
modes provided by virtual transitions. 



2 Distribution functions. Evolution equation. 

Let us consider a hypothetical system of the quark-gluon excited local thermal phase in 
QCD where a canonical operator a and its Hermitian conjugate a + occur. We formulate 
the distribution functions (DF) of produced particles (gluons and quarks) in terms of 
point-to-point equal time temperature-dependent thermal correlation functions (CF) of 
two operators 

w(k,k',t;T) = (a + (k,t)a(k',t)) = 

= Tr[a + (k,t)a(k',t)e~ Hp }/Tr(e~ H ^) . 

Here, (...) means the procedure of thermal statistical averaging; k and t are, respectively, 
momentum and time variables, e~ H/3 /Tr(e~ H ^) stands for the standard density operator in 
the equilibrium and the Hamiltonian H is given by the squared form of the annihilation a p 



and creation a+ operators for bose- and fermi- particles, H = J2 P e p a p a p (the energy e p and 
operators a p , a+ carry some index p [9], where p a = 2tt n a /L, n a = 0, ±1, ±2, V = L 3 
is the volume of the system considered); (3 is the inverse temperature of the environment, 
(3 = 1/T. The standard canonical commutation relation 

a{k,t),a + {k',t)] ± = 5 3 {k-k') 

at every time t is used as usual for Bose (-) and Fermi (+)-operators. 

The probability to find two particles, gluons or quarks, with momenta k and k' in the 
same event at the time t normalized to the single spectrum of these particles is: 

R(k, k', t) = W(k, k\ t)/[W(k, t) ■ W{k\ t)\ . 

Here, the one-particle thermal DF is defined as 

W(k,t) = (b + {k,t) b(k,t)) , 

where 

b(k,t) = a(k,t) + (f>(k,t) 

under an assumption of occurrence of the random source-function <j)(k, t) being an oper- 
ator, in general. The two-particle DF W(k, k',t) looks like 

W(k,k',t) = (b + {k,t) b + (k',t) b(k,t) b(k',t)) . 

The evolution properties of propagating particles in a randomly distributed environment 
comes from the evolution equations 

idtb(k,t)+A(k,t) = F(k,t)+P , (1) 



- i 8 t b+ (k,t) + A*(k } t)=F + (k } t) + P , (2) 

where both b and b + are the special mode operators of the quark and gluon fields [10], P 
and F(k, t) stand for the stationary external force and the random one, respectively, both 
acting from the environment. The only operator F has a zeroth value of the statistical 
average, (F) = 0. The interaction of the particles considered with the surrounding ones 
as well as providing the propagation is given by the operator A(k, t) which can be defined 
as the one closely related to the dissipation force: 

A(k,t)= K(k,t-r) b(k,r) dr . (3) 

J — oo 

An interplay of quarks and gluons with surrounding particles is embedded into the in- 
teraction complex kernel K(k,t), while the real physical transitions are provided by the 
random source operator F(k,t) ( see eqns. (P and (|)). The random evolution field op- 
erator K(k, t) in stands for the random noise and it is assumed to vary stochastically 
with a 5- like equal time correlation function [10] 

(K + {k, t) K(k', r)) = 2 (vrp) 1/2 k 5(k - k') , 
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where both the strength of the noise k and the positive constant p —>■ oo define the effect of 
the Gaussian noise on the evolution of quarks and gluons in the thermalized environment. 
The formal solution of (|l|) in the operator form in S , (3ft 4 ) {k^ — (lu — k°, kj)) is 

KK) = + <KM ' 

where the operator a(/c M ) is expressed via the Fourier transformed operator F(kfj) and the 
Fourier transformed kernel function K{k iA ) (coming from (§)) as 

a(k ll ) = F(k fl )-[K(k ll )-uj}- 1 , 

while the function ~ P ■ [Kik^) — uj]^ 1 . The random force operator F(k,t) can be 

expanded by using the Fourier integral 



F{k,t) 



+°° du) 
-oo 2tt 



-iuit 



(4) 



where the form ip^k^) ■ c(k^) is just the Fourier operator F(k^) = if} (hp) ■ c{k^) and the 
canonical operator c(fc^) obeys the commutation relation 

c(k,),c + (k' ll )] ± = 5 i (k lx -k' IM ) . 

The function ip(k^) in (|) is determined by the condition [10] 
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3 TRDD and the space-time size. 

The enhanced probability for emission of two identical particles is given by the ratio R of 
DF in 5 , (3? 4 ) as follows: 



Rb/f(k»,k';T) 



W(k.) ■ W(k> 



(5) 



where W(k^k'^T) = b(kj and W(kJ = (b + (k^) &(*„)). Using the 

Fourier solution of equation ([!]) in S(JR,^), one can get R-ratios for DF obeying to Bose- 



R b (k fl ,k^T) = l + D b (k„k'^T) 

and Fermi-particles 

R f {k^-T) = R b {k^-T) - 2 ^ 

where 



5(M • ggj 



(6) 



(7) 



D b (k^, k^j 



E(k„k'^[E(k^k,) + + 

W(k,) ■ W(k>) 



and the two-particle CF H(fc At , k' ) looks like 



- • (c + (Ag c(Ag) . (9) 



Using the Kubo-Martin-Schwinger condition ( fi is the chemical potential) 

(a(k', f) a + (k, t)) = (a + (k, t) a(k', t - i(3)) ■ exp(-/3 fi) , 

the thermal statistical averages for the c(/c /J )-operator should be presented in the following 
form: 

(c+(Ag c(k'J) = 5 4 (k,-k' fl ).n(u,T) , 

c + (^)) = 5 4 (A; M -A;;)-[l±n(^T)] 

for Bose (+)- and Fermi (-)-statistics, where n(u,T) = {exp[(o; — y)0\ ± l}^ 1 . Inserting 
CF (|^) into (H) and taking into account that the 5 4 (k^ — fc' )-function should be changed 
by the smooth sharp function Q(r) ■ exp(— q 2 /2), one can get the following expression for 
the ZVfunction 

D b (k„ T) = A(fc M , T) exp(-g 2 /2) 
x[n(u,TMr)exp(-q 2 /2) +^(k' fM )^) + fr(k^(k > fl )} , (10) 



where 



Q( r ) i 
X(k u ,k';T) = -= ^4 n(u),T) , cD = -(u; + u/) . 



The function Q(r) - n(uj; T) -exp(— q 2 /2) describes the space-time size of the QGP fire-ball. 
Choosing the z-axis along the two-heavy-ion collision axis one can put 

q 2 = (r„ • Qo) 2 + (r z • Q 2 ) 2 + (r t • Q 4 ) 2 , 

= (k — k ) jLt , Qo = e £ — e feM = &z ~ k z , Qt = — fc^,) + (fcj, — fc^) ] , 

fi(r) ~ r ■ r 2 • r 2 , 

where tq, r z and r 4 are time- like, longitudinal and transverse "size" components of the 
QGP fire-ball. Formally, the function ([10]) is the positive one ranging from to 1. 
The quantitative information (longitudinal r z and transverse r t components of the QGP 
spherical volume, the temperature T of the environment) could be extracted by fitting 
the theoretical formula ( |10"1) to the measured TRDD function and estimating the errors of 
the fit parameters. Formula (|T0| ) indicates that a chaotic g/q source emanating from the 
thermalized g/q fireball exists. Hence, the measurement of the space-time evolution of the 
g/q source would provide information of the g/q emission process and the general reaction 
mechanism. In formula ( |10| ) for the D^- function, the temperature of the environment 
enters through the two-particle CF H(fc |4 , k'^; T). If T is unstable the ^//-functions (||) will 
change due to a change of DF W which, in fact, can be considered as an effective density 
of the g/q source. Formula (|B|) looks like the following expresion for the experimental 
R-ratio using a source parametrization: 

R T (r) = 1 + Ar(r) • exp(-r 2 • Q 2 /2 - r 2 • Q 2 /2) , 



where r t (r z ) is the transverse (longitudinal) radius parameter of the source with respect 
to the beam axis, At stands for the effective intercept parameter (chaoticity parameter) 
which has a general dependence of the mean momentum of the observed particle pair. 
Here, the dependence on the source lifetime is omitted. Since < At < 1, one can 
conclude that the effective function At can be interpreted as a function of the core particles 
to all particles produced. The chaocity parameter At is the temperature-dependent and 
the positive one defined by 

= Mr) n{u-T)\ 2 
W(k») ■ W(k>J ' 
Comparing (|8|) and fllUD one can identify 

S(A^ fey = fl(r) ■ n(oo; T) ■ exp(-g 2 /2) . 

There is no a satisfactory tool to derive the precise analytical form of the random source 
function 0(fe M ) in (J|), but one can put (see @ and taking into account 0(fe M ) ~ P ■ 
[K^-uj]- 1 ) [11,10] 

0(fe,) = [«-S(fe M )] 1/2 , 
where a is of the order O (P 2 /n(u,T) ■ \ip(k^)\ 2 ). Thus, 



D b (q 2 ; T) = n f 2 g ; p. e-^ [AV^ ; T )e-^ + 2(aa'f 2 } . (11) 
(l + aj(l + a'J L J 

It is easily to see that in the vicinity of q 2 ~ one can get the full correlation if a = a' = 



and A (a;; T) = l. Putting a = a' in ( |lT|) we find the formal lower bound on the space-time 



dimensionless size of the fire-ball for Bose-system: 

A>;T) 



qt > ln 



2 



[^(q + l) 2 + a 2 - a] 
In the case of Fermi-particles, the following restriction on q 2 is valid (see (|7])) 

, Hu;T) ^2^1 A(w;T) 

ln — ^ g < Qf < m — :== 2 • 

[J2a(a + 1) + 3 - a] ' [Va 2 + 2 - a] 



In fact, the function D b (k^, k'^ T) in (0) could not be observed because of some model 
uncertainties. In the standard consideration, the TRDD-function has to contain a back- 
ground contribution as well as other physical particles (resonances) which have not been 
included in the calculation of the D b - function. In order to be close to the experimen- 
tal data, one has to expand the D b - function as projected on some well-defined function 
(in S'(3?4)) of the relative momentum of two particles produced in heavy-ion collisions 
P ) b{k IJL1 k' fl ]T) — > D b (Q 2 -T). Thus, it will be very instructive to use the polynomial ex- 
pansion which is suitable to avoid any uncertainties as well as characterize the degree of 
deviation from the Gaussian distribution, for example. In (— oo, +oo), a complete orthog- 
onal set of functions can be obtained with the help of the Hermite polynomials in the 
Hilbert space of the square integrable functions with the measure dfj,(z) = exp(— z 2 /2)dz. 
The function D b corresponds to this class if 



/+oo 
dgexp(-g 2 /2) \D b {q)\ n <oo,n = 0,l,2, 
-co 



The expansion in terms of the Hermite polynomials H n (q) 

D b {<l) = A E c n ■ H n {q) ■ exp(-g 2 /2) (12) 

n 

is well suited for the study of possible deviation from both the experimental shape and 
the exact theoretical form of the TRDD function D b (|10|). The coefficients c n in ([12]) 
are defined via the integrals over the expanded functions D b because of the orthogonality 
condition 

/+oo 
H n (x) H m (x) exp(-x /2) dx = 5 n , m . 
-oo 

Thus, the observation of the two-particle correlation (both for Bose- and Fermi-symmetrization) 
enable to extract the properties of the structure of q 2 , i.e. the space-time size of QGP 
formation. 

In order to be close to an experiment one has to replace Rj,j functions (|6]), (0) with 
respect to the cylindrical symmetry angles 9 and which are non-observable ones at fixed 

Qt- 

R bJ (k^ fcj,; T) -> R b j(Q t ; T) = N~ l J dq t dQ z d6 d<\> W(k^ T) , 

where 

N= [ dq t dQ z dO d(j> W(kp) W(k'^) , 



qt = „ 1 . „ \k x + k y T \Qt [cos(0 + 4>) + sin(0 + 4> 

cos d + sin I 2 

Then, R bJ (Q t ;T) = 1 + D bJ (Q t ;T) with 

Db ®*> T ) = 7TZ u I a exp[-(r,g,) 2 ] i*(Q t ; T) , 

D f(Qf> T ) = 7TZ via. * {#-Wexp[-(r t Q0 2 ]iW*;T) -2} , 
(1 + a)(l + a'J L J 

F b {Q t ;T) = J dq t dQ z d6d(j)n 2 {u!;T) e~ A)Z [1 + 2^ aa'A- 1 (w; T) e q2/2 } , 

/3oz = (r Qo) 2 + (r,Q,) 2 , 
jV(T) = /" rfg t c/Q 2 d# d0 n(uj] T) n(w'; T) . 
To avoid the trivial result, one should restrict the a-parameter as 

a > ~X~ 1/2 (lu; T) e^ 2 - ~^\ 1/2 (lu; T) e^ 2 . 

In conclusion of this section, the experimental data are quite desired to make the analysis 
of the particle fluctuation via TRDD. 
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4 Conclusion 



1. We investigated the finite temperature DF (of produced identical particles, gluons and 
quarks) which can be both useful and instructive to infer the shape of the gluon/quark 
source-emitter. In fact, we have presented the method of extracting the intercept and 
source parameters from the shape of the TRDD-function. 

2. The relations between the CF H(fc M , k') © and the full .R-functions for Bose (H)- and 
Fermi (0)-particles at the stage of the freeze-out are obtained. We have shown the sensi- 
tivity of the correlation functions to the space-time geometry of the source-emitter (|IT)|). 
The TRDD- function Db describes the size and shape of the space-time domain where the 
secondary observed particles are generated. 

3. One can conclude that formally, the deconfmed phase size scale can be determined 
by the evolution behavior of the field operators and the critical temperature T = T c (see 
formula (0)). 

4. Since, the TRDD-function Db is the positive one and restricted by 1, we expect that the 
i?-ratio at too small values of starts from the fixed point R{Q^ — > 0) = 2 — e (e — > +0) 
and then falls down (with the Gaussian shape) up to unity over some momentum scale 
interval of an order of the inverse source size. 
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